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Abstract: Using the 3-year observational data from ChinaFlux (Chinese 
Terrestrial Ecosystem Flux Research Network), we studied the gas regu¬ 
lation flux dynamics and cumulative process of gas regulation value in 
Qianyanzhou middle subtropical plantation (QYF) and Changbai Moun¬ 
tain temperate mixed forest (CBF). The gas regulation service was dif¬ 
ferentiated into vegetation gas regulation service and net ecosystem gas 
regulation service. Carbon tax approach, reforestation cost approach and 
industrial oxygen approach were employed to calculate gas regulation 
value. Results show that there was significant seasonal variation in vege¬ 
tation gas regulation flux. Daily C0 2 uptake fluxes averaged 82.00 
kg-ha'^d' 1 and 59.37 kg-ha'^d' 1 and the corresponding 0 2 emission fluxes 
were 59.65 kg-ha'^d' 1 and 43.19 kg-ha'^d' 1 for QYF and CBF, respec¬ 
tively. The cumulative curves of vegetation gas regulation value always 
followed a sigmoid shape, and the annual gas regulation value produced 
by vegetation was RMB 14 342.69 yuan-ha’ 1 and RMB 10 384.18 
yuan-ha’ 1 for both QYF and CBF, respectively. In terms of monthly net 
ecosystem gas regulation service, QYF appeared as a C0 2 sink and 0 2 
source for the whole year, while CBF appeared to be a C0 2 sink and 0 2 
source mainly in the period between May and September. The cumula¬ 
tive curves of net ecosystem gas regulation value presented a sigmoid 
("S") shape for QYF, while a unimodal type curve for CBF. The annual 
net ecosystem gas regulation value was 8470.52 yuan-ha' 1 and 5091.98 
yuan-ha' 1 for QYF and CBF, respectively. The economic value of both 
the vegetation gas regulation service and net ecosystem gas regulation 
service were mainly produced between May and October. 
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Introduction 

As the main body of global terrestrial ecosystem, forest ecosys¬ 
tem produces important ecosystem services (Costanza et al. 1997; 
Myers 1997). In forest ecosystem, plants transform solar energy 
into biotic energy through photosynthesis, fixing C0 2 and re¬ 
leasing 0 2 , which plays an irreplaceable role in maintaining the 
C0 2 /0 2 balance and mitigating greenhouse gases emissions. 
Annually, two thirds of the global terrestrial carbon sequestration 
is achieved by forest ecosystem (Kramer 1981). Costanza et al. 
(1997) grouped the ecosystem services into 17 major categories 
for 16 biomes, in which gas regulation service is the most im¬ 
portant one. According to the equivalent weight factor of eco¬ 
system services per hectare of forest ecosystem in China, gas 
regulation service was given the high equivalent factor of 3.5, 
which is the second largest ecosystem service delivered by forest 
(Xie et al. 2003). The economic value of carbon fixation and 0 2 
emission was estimated at RMB 1626.76><10 8 yuan-ha" 1 -a" 1 and 
6732.48><10 8 yuan-ha" 1 -a" 1 for China forest ecosystem (Zhao et al. 
2004). A number of other studies also evaluated gas regulation 
services at the country and regional level (Ouyang et al. 2004; 
Jin et al. 2005; Yu et al. 2005a; Yu et al. 2005b). 

However, to date, studies on gas regulation service of forest 
were mainly based on biomass and productivity survey, or forest 
inventory (Zhao et al. 2004; Jin et al. 2005; Yu et al. 2005b). 
Hence, these studies on gas regulation service were confined to 
evaluating and analyzing variation of gas regulation service 
among different forest ecosystems at the annual or interannual 
scale, which resulted in deficiency of the revelation of seasonal 
dynamics of gas regulation service. In fact, forest ecosystem is a 
dynamic system, and the gas regulation service also presents 
seasonal variation. Therefore, a study on seasonal dynamics of 
gas regulation service contributes to revealing the formation 
mechanism of gas regulation service and its process. Moreover, 
what the previous study evaluated is virtually vegetation gas 
regulation service (Zhao et al. 2004; Jin et al. 2005; Yu et al. 
2005b), which overlooked heterotrophic respiration and resulted 
in overestimation of ecosystem gas regulation service . 

In this paper, we differentiated the gas regulation service into 
vegetation gas regulation service and net ecosystem gas regula- 
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tion service, and the former was computed based on NPP (net 
primary productivity), while the latter was calculated based on 
NEP (net ecosystem productivity). Eddy covariance system can 
continuously and directly measure the C0 2 flux between vegeta¬ 
tion and atmosphere, which meet the demand for study of sea¬ 
sonal dynamics of these two gas regulation services. The ques¬ 
tions addressed here were: (1) characteristics of the gas regula¬ 
tion flux dynamics; (2) cumulative process of economic value of 
gas regulation service. 

Materials and methods 

Site description 

The flux measurements were made above the plantations in 
Qianyanzhou Experimental Station (QYF) and broad-leaved 


mixed forest in Research Station of Changbai Mountain Forest 
Ecosystems (CBF). General situation of these two sites was 
listed in Table 1. More detailed descriptions of the two sites are 
provided by Zhang et al. (2006a), Liu et al. (2006) and Zhang et 
al. (2006b). The equipment in these two stations is uniform 
open-path eddy covariance system, and the description of the 
flux observation tower can be referred to the ChinaFlux website: 
http://www.chinaflux.org . 

Data acquisition and processing 

NEE (Net ecosystem carbon exchange), Re (Ecosystem respira¬ 
tion) and GEE (Gross ecosystem carbon exchange) flux at daily 
scale were provided by Chinese Terrestrial Ecosystem Flux Ob¬ 
servational Research Network (ChinaFlux). 


Table 1. General situation of the study sites 


Site 

Geographic 

location 

Climate 

Annual mean 

temperature (°C) 

Annual pre¬ 

cipitation (mm) 

Crown 

height (m) 

Dominant species 

Soil type 

Observational 

height (m) 

QYF 

26°44'N, 

Subtropical mon- 

17.9 

1542.4 

11 

Finns massonian, Pinus elliottii, 

Red soil 

39 

115°03'E 

soon climate 




Cunninghamia Lanceolata 



CBF 

41°29'N, 

Temperate conti- 

3.6 

713 

26 

Pinus koraiensis, Tilia turn, Quercus 

Dark brown 

40 

128°05'N 

nental climate 




mongolicus, Fraxinus mandshurica 

forest soil 



Notes: QYF—Qianyanzhou middle subtropical plantation, CBF—Changbai Mountain temperate mixed forest. 


The measurements were made in 2005, 2006 and 2007. Net 
ecosystem carbon exchange between forest and atmosphere 
(F NEE ) was computed according to the following equation: 


F 


NEE 



d c 
~dt 


dz + (w'c')r 



where, the first term on the equation is the storage term (F s ), the 
second term is the eddy C0 2 eddy flux ( F c ) measured at height Z r . 
Subscript r denotes a quantity at eddy flux measurement height 
Z r . F nee is decomposed into F Re (flux of ecosystem respiration) 
and F gee (flux of gross ecosystem carbon exchange) according to 
the following equation: 


= F™ - K 


GEE 


NEE 


Re 



More detailed data processing are referred to Zhang et al. 
(2006a), Zhang et al. (2006b). The non-linear regression tech¬ 
nique is used as gap-filling method. 

Based on flux data, C0 2 uptake and 0 2 emission fluxes were 
computed according to the formula of photosynthesis and respi¬ 
ration (Guo et al. 2001; Xiao et al. 2005b): 

6/iC0 2 +6wH 2 0^/iC 6 EIi 2 0 6 +6w0 2 —>^C 6 H 10 O 5 

When 1-g carbon was sequestered, the 3.67-g C0 2 is fixed and 
2.67-g 0 2 is released. On ecosystem scale, GPP (Gross primary 
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productivity) is the modulus of GEE. And GPP minus plant res¬ 
piration equals NPP. Assuming the ratio of NPP/GPP (the carbon 
use efficiency) kept unchanged in the year, the flux of C0 2 up¬ 
take and 0 2 emission by vegetation could be determined by the 


equation as follows: 


A c = “Nee x / x 3 .67 x 10 

(3) 

A o = -Nee x / X 2.67 X 10 

(4) 


where, A c is the C0 2 uptake flux by vegetation (kg-ha'^d 1 ), A 0 is 
the 0 2 emission flux by vegetation (kg-ha^d 1 ), F GEE is the flux 
of gross ecosystem carbon exchange (gC-m^-d 1 ),/^ the carbon 
use efficiency, and 10 is a conversion factor for gC-m' 2 to kg-ha' 1 . 
The carbon use efficiency for QYF and CBF was determined by 
referring to literatures (Ma et al. 2008; Wang et al. 2006; Zhang 
et al. 2008). 

NPP quantifies the carbon uptake only by plants, while NEP 
(negative NEE) includes carbon absorption by plants and carbon 
release by soils. Properly speaking, NEP is net carbon exchange 
between the ecosystem and the atmosphere. Therefore, the gas 
regulation service by ecosystem should be determined by NPP 
rather than NPP, and it was called net ecosystem gas regulation 
service in this paper. Net ecosystem gas regulation service was 
defined as the following equations: 

N c =-F nee x3.67 x 10 (5) 
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Nq —F NEE x 2.67 ><10 (6) 

where, N c is the net ecosystem C0 2 uptake flux (kg-ha^-d' 1 ), N 0 
is the net ecosystem 0 2 emission flux (kg-ha''-d''), F NEE is the 
flux of net ecosystem carbon exchange (gC-m^-d' 1 ), and 10 is a 
conversion factor for gC-m' to kg-ha' . 

The economic value of C0 2 uptake was estimated by the av¬ 
erage value of the cost of afforestation and Swedish carbon tax 
(Ouyang et al. 1999; Xiao et al. 2005a), i.e. RMB 0.20535 
yuan-kg' 1 C0 2 . The economic value of 0 2 emission was evaluated 
using the average value of the price of industrial 0 2 (0.4 
yuan-kg" 1 ) and the cost of afforestation (RMB 0.3529 yuan-kg' 1 ) 
in China (Ouyang et al. 1999; Xiao et al. 2005a), i.e. 0.37645 
yuan-kg' 1 0 2 . 

The cumulative economic value of gas regulation service was 
calculated as follows: 

A, = t v, (7) 

-1 

where, represents the cumulative value of gas regulation ser¬ 
vice for i months, and V t represents the economic value of gas 
regulation service in zth month. In the paper, the economic value 
of C0 2 uptake and 0 2 emission was calculated separately. 


Results and discussion 

Vegetation gas regulation service 

Over the 3 years, the NPP were (815.50±55.18) gC-m^-a" 1 and 
(590.43±25.50) gC-m^-a' 1 for QYF and CBF, respectively. These 
results agree well with the earlier studies (Ma et al. 2008; Wang 
et al. 2006; Zhang et al. 2008). Ma et al. (2008) found that NPP 
simulated by BGC model ranged from 343.31 gC-m^-a" 1 to 
906.42 gC-m'^a' 1 across twenty-one years for QYF. As to CBF, 
the NPP value was also consistent with results of Zhang et al. 
(2008) and Wang et al. (2006). Due to the seasonal variation of 
NPP, vegetation gas regulation service showed obvious seasonal 
dynamics. The seasonal variation of vegetation gas regulation 
service for both QYF and CBF were unimodal distribution (Fig. 
1). The 3-year average daily C0 2 flux of QYF and CBF is 
(82.00±5.55) kg-ha^-d" 1 and (59.37±2.57) kg-hm'^d" 1 , respec¬ 
tively. Correspondingly, the daily average 0 2 flux was 
(59.65±4.04) kg-hm^-d^and (43.19±1.86) kg-ha'^d" 1 , respec¬ 
tively. From January to April, C0 2 fluxes of both QYF and QBF 
were lower than their daily average C0 2 fluxes, which should be 
attributed to the lower temperature in environment and lower 
growth rate of vegetation. The peak value of C0 2 uptake flux and 
0 2 emission flux of QYF occurred in August, while CBF in July, 
which may be ascribed to the climate characteristics difference. 



1 2 3 4 5 6 7 8 9 10 11 12 

Month 



Fig. 1 Flux of C0 2 uptake (a) and 0 2 emission (b) by vegetation 

CBF— Changbai Mountain temperate mixed forest; QYF— Qianyanzhou middle subtropical plantation. 


With the biomass formation, the economic value of vegetation 
gas regulation service is produced and accumulated. The cumu¬ 
lative dynamic curves of vegetation gas regulation value showed 
a sigmoid shape for QYF and CBF, and the annual gas regulation 
value produced by vegetation was (14342.69±970.54) yuan-ha' 1 
and (10384.18±448.51) yuan-ha' 1 , respectively (Fig. 2). In our 
study, annual economic value of C0 2 uptake averaged 6145.91 
yuan-ha' 1 and 4449.67 yuan-ha' 1 for QYF and CBF, respectively. 
These results are similar to that of Ouyang et al.(2004), who 
found that the economic value of C0 2 uptake for tropical mon¬ 
soon forest was 5343.3 yuan-ha' 1 . The value accumulated in 
spring, summer, autumn and winter accounted for 16.83%, 
36.44%, 33.39% and 13.34% of the annual total value for QYF, 


respectively. However, the value accumulated in summer and 
autumn occupied 52.07% and 38.07% of the annual total value 
for CBF, respectively. Additionally, the vegetation gas regulation 
service mainly occurred between May and October, and the 
economic value produced during this period took up 69.83% and 
90.14% of the annual total value for QYF and CBF, respectively. 
Besides, the monthly vegetation gas regulation value of QYF 
was always larger than that of CBF. 

Net ecosystem gas regulation service 

For the whole forest ecosystem, net ecosystem gas regulation 
service equals to vegetation gas regulation service subtracted by 
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heterotrophic respiration. The 3-year average NEP was 
(481.62±54.38) gOm'V 1 and (289.52±41.97) gC-m’V 1 for 
QYF and CBF, respectively, and the results are consistent with 
Zhang et al. (2006a). Moreover, there were significant seasonal 
trend in the fluxes of C0 2 uptake and 0 2 emission for both QYF 
and CBF (Fig. 3). At the monthly scale, QYF served as C0 2 sink 
and 0 2 source for a whole year, and similar results were also 
reported in other studies (Liu et al. 2005; Zhang et al. 2006a). 
CBF served as a C0 2 sink and 0 2 source mainly in the period 
between May and September, while in other months CBF main¬ 
tained rough C0 2 balance or acted as a weak C0 2 source. Daily 
average C0 2 fluxes for QYF and CBF were 48.43 kg-ha'^d' 1 and 



Month 


29.11 kg-ha'^d" 1 , respectively, and daily average 0 2 fluxes were 
35.23 kg-ha'^d' 1 and 21.18 kg-ha'^d' 1 for QYF and CBF, respec¬ 
tively. The largest net ecosystem C0 2 uptake flux and 0 2 emis¬ 
sion flux for both QYF and CBF occurred in July, and the largest 
net ecosystem gas regulation flux of CBF was 1.73 times as large 
as that of QYF. Liu et al. (2005) observed that the highest NEE 
(negative NEP) occurred in May or June, which differ from our 
results. It was attributed to severe drought in QYF in 2003. Gu et 
al. (2008) also demonstrated that drought affect the NEP. Addi¬ 
tionally, the seasonal variation of C0 2 uptake flux and 0 2 emis¬ 
sion flux is greater than that of CBF. 



Month 


Fig. 2 Cumulative process of economic value of CO 2 uptake (a) and 0 2 emission (b) by vegetation 



Fig. 3 Net ecosystem C0 2 uptake (a) and 0 2 emission (b) fluxes 




Month 


Fig. 4 Cumulative process of economic value of net C0 2 uptake (a) and net 0 2 emission (b) by ecosystem 


There was significant difference in cumulative process of QYF and CBF (Fig. 4). The cumulative curves of net ecosystem 

economic value of net ecosystem gas regulation between the gas regulation service for QYF presented sigmoid type curve, 
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while that for QYF followed a unimodal shape. The annual eco¬ 
nomic value of C0 2 uptake and 0 2 emission was 8470.52 
yuan-ha" 1 for QYF, which was 1.66 times as large as that of CBF. 
Xiao et al. (2005b) reported that the value of 0 2 emission for 
subtropical rice paddy ecosystem ranged from 9549 to 12277 
yuan-ha" 1 , which is much larger than that of QYF and CBF, and 
we supposed that it was attributed to higher productivity of rice 
paddy ecosystem. Between May and October, the cumulative 
economic value of net ecosystem gas regulation of QYF and 
CBF accounted for 66.96% and 108.70% of the annual total 
value. The economic value of net ecosystem gas regulation for 
CBF reached maximum in September, 5602.91 yuan- ha" 1 , but in 
October and November CBF became a C0 2 source and 0 2 sink, 
which brought about negative economic value and decreased the 
cumulative economic value. 

Conclusions 

There were significant seasonal variation of vegetation gas regu¬ 
lation fluxes for QYF and CBF, and the dynamic curves of 
vegetation gas regulation fluxes for QYF and CBF presented 
unimodal shapes. Annual economic value of vegetation gas 
regulation for QYF and CBF was RMB 14342.69 yuan-ha" 1 and 
10384.18 yuan-ha" 1 , respectively, and their cumulative dynamics 
followed sigmoid type curve. The seasonal variation of net eco¬ 
system gas regulation service fluxes for QYF and CBF was also 
significant. In terms of monthly net ecosystem gas regulation 
service, QYF served as a C0 2 sink and 0 2 source for a whole 
year, while CBF served as a C0 2 sink and 0 2 source for eight 
months. The cumulative dynamics of economic value of net 
ecosystem gas regulation appeared sigmoid curve for QYF and 
unimodal curve for CBF. Annual economic value of net ecosys¬ 
tem gas regulation was 8470.52 yuan-ha" 1 and 5091.98 yuan-ha" 1 , 
respectively. 
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